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A theoretical description of scanning tunneling potentoimetry (STP) measurement is presented to 
address the increasing need for a basis to interpret experiments on macrscopic samples. Based on a 
heuristic understanding of STP provided to facilitate theoretical understanding, the total tunneling 
current related to the density matrix of the sample is derived within the general framework of 
quantum transport. The measured potentiometric voltage is determined implicitly as the voltage 
necessary to null the tunneling current. Explicit expressions of measured voltages are presented 
under certain assumptions, and limiting cases are discussed to connect to previous results. The 
need to go forward and formulate the theory in terms of a local density matrix is also discussed. 
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I. INTRODUCTION 

The use of scanning probes is now widespread. How- 
ever, very few of these probes measure transport at very 
short length scales. The ultimate probe for such measure- 
ments is scanning tunneling potentiometry (STP)^ ^, in 
which a scanning tunneling microscope (STM) is used 
to measure the local potential due the flow of an ap- 
plied current. However, only recently have we and oth- 
ers developed STP instruments that can routinely make 
measurements under a wide range of conditions and es- 
sentially at the fundamental noise limit of STP.E1E1_ 

With our STP, we have focussed on local quantum 
transport in macroscopic materials, as distinct from 
quantum transport through nanostructures (e.g., sin- 
gle molecules, nanotubes and lithographically produced 
nanostructures). Our work demonstrates that potential 
maps can be obtained in appropriate materials down to 
distances smaller than all the length scales relevant in 
transport (the inelastic scattering length, the elastic scat- 
tering length and even the Fermi wave length). An ex- 
ample of an STP image on these very short length scales 
obtained in our laboratory is shown in Figure [T] for an 
epitaxial graphene sample. As can be seen in the figure, 
there is considerable local structure in the potential that 
is relatively large in magnitude. It is the lack of a quan- 
titative basis for interpreting such images that motivates 
this paper. We hope both to explain the theoretical issues 
for experimentalists and to motivate the needed theoret- 
ical extensions for theorists. 

There are two issues. First, it is not clear exactly what 
potential is being measured. Macroscopically, the mea- 
sured potential would be the usual local electrochemical 
potential. However, at such short length scales as STP 
probes, one cannot use thermodynamic concepts. Sec- 
ond, the situation is inherently quantum mechanical and 
any calculation of the potential must include the under- 
lying quantum mechanical processes. 

These issues are not entirely new. To some degree they 
have been addressed in the theories relevant to trans- 
port through nanostructures where the sample is small 



compared to characteristic lengths of the transport. On 
the other hand, in the case of these nanostructures, in 
most of the cases the measurement contacts are relatively 
macroscopic (compared with an STP tip), and cannot 
be scanned. Moreover, the interfaces between the con- 
tacts and the sample play an important role in the over- 
all transport. Numerous theories have been developed 
to deal with this situation. Some of them also consider 
what an STP would measure inside a nanostructure2"m, 
although there are no existing STP data yet that can be 
compared with these theories. In any event, regarding 
the specific question of what STP measures, it was possi- 
ble to write down explicit expressions for the STP voltage 
(for example, see equation (37) of reference [7]) . On the 
other hand, these expressions are explicitly dependent on 
the voltages applied on the current leads. By contrast, in 
the case where the sample is macroscopically large, the 
geometry and microscopic processes present in the leads 
obviously cannot matter. Hence a proper theoretical for- 
mulation is required, and precisely what STP measures 
remains unclear. 

In this paper, we take some first steps toward a theoret- 
ical description of a STP measurement of local quantum 
transport in a macroscopic sample. The relevant quan- 
tum mechanical quantity is the density matrix, and it 
is possible to develop an implicit (and under some ap- 
proximations, an explicit) relation from which the mea- 
sured potential can be determined in terms of the den- 
sity matrix and the tunneling process into the sample. 
While the tunneling process associated with the STM 
tip can be accounted for, it is a separate matter how to 
do the calculation of the density matrix under the rel- 
evant nonequilibrium conditions and short length scales 
in a macroscopic sample. We do not address these com- 
putational challenges in this paper. As stated above, our 
goal is to understand scanning tunneling potentiometry 
in terms accessible to experimentalists and to define the 
deeper theoretical questions needing attention. Toward 
these ends, we use both heuristic and more formal ap- 
proaches. More specifically, in the formal development, 
our results pose new problems in quantum transport that 
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FIG. 1. STP data taken on epitaxial graphene as an example 
of transport measurement at length scales smaller than all the 
length scales relevant in the transport. As seen in the figure, 
on top of the average gradient, STP shows variations on the 
nanometer length scale. 



need proper theoretical attention in order to have a com- 
plete theory suitable for interpreting experiment. One 
example is the need for formulations in terms of local 
density matrices as opposed to global ones. 

We note that some of the issues discussed above have 
been addressed by Chu and Sorbellc^^ in the context 
of their calculation of the residual resistivity dipole of 
a scattering center in the spirit of Landauei^^. As noted 
in their paper, the potential measured in STP is not 
the same as the local electrostatic potential in the mate- 
rial generated by an electric current. Our work can be 
thought of as a generalization of their work and an artic- 
ulation more from the general point of view of a theory 
of STP measurement. Compared to Chu and Sorbello, 
our result is more general; and in the final result, instead 
of a weighted sum, we obtained a form which consists 
of matrices either being density matrix of the sample or 
defined with sample and tip wave functions. 

This paper is organized as follows. First, we describe 
the setup for STP measurement, also serving to empha- 
size that what is operationally measured is the potential 
necessary to achieve zero current through the potentio- 
metric contact. We provide a heuristic understanding of 
STP measurement in general terms. Next, we present a 
simple model of only two scattering centers with meso- 
scopic distance to motivate a more general theory. This 
leads naturally to the use of the density matrix of the 
sample in the expression of the tunneling current de- 
tected in STP and thereby gives meaning to what is being 
measured. After that, we discuss the need to use a local 
density matrix in the case of macroscopic samples. Fi- 
nally, we discuss several limiting cases of our theory to 
illustrate better its physical content. 




FIG. 2. Schematic setup of scanning tunneling potentiometry, 
and model problem to be considered in this paper. Electrodes 
1, 2 and 3 are macroscopic contacts, while the STM tip, mi- 
croscopically connected to the sample, functions as the fourth 
electrode. The sample is macroscopic in size and STP probes 
a small region which is enlarged to the upper left. The re- 
gion has two scattering centers, and the rest of the sample is 
assumed to be homogeneous and defect-free, with no defects 
(scatterers) present. The distance between the two defects, 
the Fermi wavelength and the inelastic mean free path of the 
sample are all comparable. To the upper right presents an 
expected STP result on a large length scale, including the 
electrodes. 



II. SCANNING TUNNELING 
POTENTIOMETRY 



A. STP setup 



As seen in Figure [2j STP is effectively a four-point 
transport measurement using an STM tip. In an STP 
measurement, a floating current source provides a cur- 
rent through the sample via electrodes 1 and 2; a volt- 
age is applied between the third electrode and the STM 
tip, which serves as the fourth electrode. The voltage 
is so adjusted that the tunneling current between the 
sample and the STM tip is zero - the definition of a po- 
tentiometric measurement. This applied voltage which 
nulls the tunneling current is the data that is recorded in 
an STP measurement. Moreover, the capability of STM 
to scan on nanometer scales makes STP measurement a 
nanoscale transport measurement. The outcome of the 
experiment is in the form of a potential map of the scan- 
ning area accompanied by a topographical map (obtained 
by conventional STM operation) of the same area taken 
point by point successively with the potential (see Ref |3I 
for details). 



B. A heuristic interpretation of STP 

Now consider a macrocopic sample where the current 
contacts are far removed, specifially where the sample 
size and the distance between the current contacts are 
large compared to characteristic transport length scales. 
Our goal is to address the question what STP really mea- 
sures under these conditions. This measured value in 
STP would correspond to the local electrochemical po- 
tential in the thermodynamic limit (i.e., a large distance 
between the voltage measurements and a large voltage 
contact area in addition to a macroscopic sample size). 
In STP the situation is more subtle. 

In a conventional STM measurement, the total tunnel- 
ing current is writterP^ 
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where 5{Ei,E2) is the Kronecker delta symbol, fi (i — 
s, t) are Fermi-Dirac distribution functions of the sample 
and the tip, respectively; |Mp is the magnitude of tun- 
neling matrix elements, which are assumed to be contant 
in this equation; E^^ g are eigen energies of the specific 
states of the sample and the tip, where the energies here 
are measured with respect to the band ("chemical en- 
ergy" , see Figure |3] for clarifying definitions of chemical 
energy, electrostatic energy and the sum of these two, the 
electrochemical energy); and —eVt are the energies asso- 
ciated to the electrostatic potential in the sample and the 
tip ("electrostatic energy", also see Figure [s]). 

Two modifications need to be made in order to use 
the form of equation (IT]) in the case of a potentiometric 
transport measurement. First, the distribution function 
of the sample is not one in equilibrium due to existence of 
an applied current through the sample, and hence is not a 
Fermi-Dirac distribution. And second, the magnitude of 
tunneling matrix elements |Mp should not be assumed 
as a constant for each tunneling channel, hence in the 
expression one should use an equivalent average value 
for the magnitude of tunneling matrix elements. 

Consider now an example of the distribution function 
in an non-equilibrium transport situation. In the linear 
response region of a homogeneous sample with no defects 



and only inelastic scattering with mean free path /i„, the 
set of electrons with a certain direction of wave vector 
is described by an effective chemical potential dependent 
on this directioiPl, i.e., fs{k) ex {el^(^-t's) + l)-i (here 
(3 = 1/(/cbT), see also Section |v]). A distribution func- 
tion fs{Es^ec) can be defined as the average probability 
of occupancy for each energy Es^ec, where the averaging 
runs over all degenerate states (with different direction of 
wave vector, for example). Plotting fs{Es.ec) as a func- 
tion of Es_ec, it is expected that /s ~ 1 at low energy 
and /s ~ at high energy, but in the transition region 
between low and high energies, it is not a Fermi-Dirac 
distribution. If we ignore the broadening of the distribu- 




FIG. 3. Definitions of different energies in the problem. 
Chemical energies (Ec) are defined with respect to the band, 
e.g., the band bottom; electrostatic energies (Ee) are energy 
differences associated with the band bending due to existence 
of an electrostatic field; and the sum of these two energies 
are the electrochemical energies (-Eec). Electrons at different 
positions of the sample with the same chemical energy and 
wave vector have the same wave function, while they do not 
necessarily have the same electrochemical energy due to band 
bending. Although electrochemical energy is the appropriate 
energy one needs to use when comparing energies at different 
positions, it is easier to keep track of quantum states with 
chemical energies. 



tion function due to finite temperature, with a constant 
current and a cylindrical Fermi surface, in the previous 
case where the sample is homogeneous and defect-free, 
the distribution function is in the form (see also Figure 
HI for more detailed clarification): 
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1 ^s,ec < M ~ ^hnEo 



fs{Es,ec) — < 
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2{ehnEof/^ 




fl - eknEo < Es^ec < fJ- 

M < Es^ec < /^ + elinEo 
Es,ec > M + elinEo 
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FIG. 4. Comparison of distribution functions with and with- 
out current, (a) When there is a current, the equihbrium 
Fermi surface (dashed hne and shaded) is shifted to the right 
(solid Hne). (b) When viewed in the E — k diagram, in equihb- 
rium, the Fermi surface is a horizontal circle, (c) while when 
current exists, the Fermi surface is tilted, (d) The occupa- 
tion rate for each energy, i.e., the distribution function (solid 
red curve), is different than Fermi-Dirac distribution (blue 
curve). The distribution function in more complicated situa- 
tions, for example, when elastic scattering exists, is expected 
to be different and more interesting than illustrated here. 



nitude of the tunneling matrix represents the quantum 
interference, which we will study in detail in next sec- 
tion. This concept that the total current is proportional 
to a weighted sum of the difference of two distribution 
functions is not unlike thermal electric effect, although 
in that case it is the difference in temperature that is the 
origin of the difference in distribution function, while in 
our case it is the nonequilibrium nature of current flow. 
From equation (Is]) it is clear that STP does not mea- 
sure a well defined thermodynamic potential, rather, it 
measures the non-equilibrium distribution function of the 
sample with the weight for each energy due to both den- 
sities of states and quantum interference effects in the 
tunneling matrix. To actually calculate this weight, a 
density matrix is needed, which is described below, ft is 
also clear that even with the same distribution function, 
due to relative change in the tunneling matrix element 
from quantum interference from one measurment point 
to another, the STP voltage can change. This accounts 
for the STP fluctuation in Chu & Sorbello's papeM The 
basic message of this heuristic consideration is that the 
STP potential reflects both changes in the distribution 
functions and the tunneling matrix elements as a function 
of position. These insights provide a qualitative basis for 
considering experiment and set the stage for a more exact 
theoretical treatment. 



III. GENERAL THEORY IN TERMS OF 
DENSITY MATRIX 

A. Simple model considered in the problem 



where E^ is the electric field and lin is the inelastic mean 
free path. 

Comparing the distribution function of the sample to 
that of the tip (a Fermi-Dirac distribution), it is clear 
that electrons with relatively low energies will tend to 
flow from the tip to the sample, and electrons with rela- 
tively high energies will tend to flow in the opposite direc- 
tion. Hence the total tunneling current is a weighted sum 
of the difference of the two distribution functions. The 
weight for each energy is clearly related to the density 
of states, however, this weight is also related to the tun- 
neling matrix elements as they are determined by wave 
function values at a certain position (tip center of cur- 
vature) and should not be assumed to be the same as in 
equation ^\. One can write: 
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(3) 
where V^, which is the voltage on the tip, effectively 
moves the distribution function ]t (blue curve in Figure 
llld)) in the horizontal direction, such that the total tun- 
neling current equals zero. Here, the densities of states 
are more or less semi-classical while the average mag- 



We introduce a simple model presented in Figure [2] 
to represent a macroscopic sample with defects closer to 
one another than the inelastic mean free path. In the 
problem, the sample is macroscopic and connected to the 
current providing electrodes macroscopically. The STM 
tip probes a small region on the sample, and inside this 
small region are two scattering centers near each other. 
The rest of the sample is assumed to be defect-free and 
homogeneous, with inelastic mean free path lin and Fermi 
wavelength \p. These two parameters are both com- 
parable to the distance between the two defects in the 
problem. The electrochemical potential on the sample 
is expected to be a linear curve with respect to position 
except for where the defects are, as also shown in the 
figure. However, as one expects from physical consid- 
eration, and explained in detail below, when comparing 
STP data taken at nearby points, such thermodynamic 
concepts are inadequate. 

In this problem, one can ask the following two ques- 
tions: a) what is the STP measurement result near the 
two defects? b) what is the cross section of the two de- 
fects, if seen far away and treating them as one single 
defect? The answer to the second question determines 
STP measurement result far away from the defects, as 
in that case STP measures the Landauer resistive dipole 



potentiapSI. One needs to solve a similar problem to an- 
swer both questions. 

B. Quantum transport approach 

In the quantum transport approaclpSl^ one uses a den- 
sity matrixes to describe the sample, and treat all three 
contacts as electron reservoirs that have Fcrmi-Dirac dis- 
tributions. The necessity of density matrix in solving the 
model problem posed above can be appreciated with the 
following consideration, with plane waves as the basis set 
for electron wave functions: 

Suppose we start with an incident plane wave. The 
plane wave will be scattered by either of the two de- 
fects. The scattered wave is coherent with the incident 
wave. As it propagates, the magnitude of the coherent 
part shrinks due to inelastic scattering, which generates 
waves having indefinite phase difference with the incident 
and scattered waves. Moreover, the scattered coherent 
wave can be scattered by the other defect, generating a 
coherent secondary scattered wave. In short, the inci- 
dent plane wave will be scattered back and forth by both 
scattering centers, with decaying coherent part. In addi- 
tion, the incoherent waves generated by inelastic scatter- 
ing will individually also be scattered back and forth and 
generate coherent sets of their own. The situation which 
combines coherence from elastic scattering and incoher- 
ence and probability from inelastic scattering is formally 
described by a density matrix. The density matrices re- 
sulting from all the incident plane waves which are inco- 
herent to one another should be added to obtain the net 
density matrix. 

In this paper, we do not undertake the calculation of 
the density matrix of the sample. Rather, we assume 
that the density matrix has been calculated. We did 
include the macroscopic contacts in the problem, since 
they are needed in determining the density matrix of the 
sample in the Non Equilibrium Green Function-Landauer 
approach^^. The objective of this paper is to calculate 
the total tunneling current based on the density matrix of 
the sample and properties of the STM tip. Once the total 
tunneling current is obtained with the voltage applied on 
the tip as a parameter, by setting that current to zero 
one implicitly determines STP measurement result. 

C. Total tunneling current 

In calculating the total tunneling current, we have the 
density matrix of the sample, denoted by g, and describe 
the STM tip with a Fermi-Dirac distribution function: 



f{E,,j;-eVt-^^) 



1 



To derive the expression for the total tunneling current, 
we start by diagonalizing the density matrix g of the 
sample . After diagonalization, the density matrix reads: 



/•. 



P-0 = 
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Pipkk 



P-Lpk'k' 



V- 



(4) 
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where P^^tj: = (V'jrl^lV'fc)- Here we use g for the abstract 
density matrix, and use p^ for the matrix form of the 
density matrix under the basis set {jf/^j:)}. In this partic- 
ular case, this basis set diagonalizes the density matrix. 

The basis set (I'/'j:)} consists of wave functions that 
are incoherent to each other in the problem, and the di- 
agonal elements in the density matrix p^ have physical 
meanings that the probability of finding an electron in 
state IV'g) is P^^^- Note that elastic scattering only scat- 
ters states with a given eigen electrochemical energy into 
states with the same eigen electrochemical energy, thus 
one can choose {|^^:)} such that they are all energy eigen 
states, with E^^ ^r as their eigen electrochemical energies. 
Since we are using chemical energies for states in the tip 
and electrochemical energies for states in the sample, the 
subscripts "s" and "t" in the denotation of energy levels 
for sample and tip have been omitted. 

As shown in Figure [5] the tunneling current between 
one state jV'j:) of the sample and one state \xp) of the tip 
is calculated through evaluating a surface integral over S 
between the sample and the tip, and reads--: 



^k,p ^ fc [Pipkk ' 
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fiE,^p-eVt 



t^MEj:,E,,j,~ 
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eVt) 



(5) 



The tunneling current resulting from different channels 
are additive, because both the states of the sample and 
the states of the tip are incoherent with one another. 
Thus the total tunneling current is: 

^^T.\^ [P^kk-fiEc.^p-^Vt~t.)]diE^^^^„E,,p-eVt) 

kn '^ 



2m 



i^k^xp - x*jNi^k) ■ d5 



(6) 



The tunneling matrix element: 



with llxp)} as its basis set. Here Vt is the voltage applied 
on the STM tip, Ef,^p is the chemical eigen energy of state 
IXp), and /^ is the electrochemical potential of the tip. 



Mi,^5{E^^^^,,E,^p^eVt)— / (^,,Vx;~ - x;-V^^) ' d^ 

(7) 



can be written in another way if one expands the wave 
function of the tip into spherical harmonics (see chapter 
3.2 of reference [H]) : 



Xp(^ = ^Ci.m,pkl{Kp)Yirn{9,(j)) 



(8) 



l,m 



where besides the famihar denotation of Yim as the spher- 
ical harmonic functions, ki is the l^^ spherical modified 
Bessel function of the second kind. There exists a gen- 
eral derivative rule (chapter 3.4 of reference [H)) that the 
density matrix element ([7| can be written as a linear op- 
eration of the sample state only, evaluated at the center 
of curvature rg of the tip which is also the origin for ex- 
panding the tip wave functions: 



%p = i^ipyt)i'k) \ro 



(9) 



For instance, if the tip wave function is well approximated 
by a s-wave, then: 
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The total tunneling current can be expressed: 
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(10) 

where the Kronecker delta symbol has been absorbed in 
the operator J^(p,Vt)j ^^'^ these operators are well defined 
only between energy eigen states. From the point of view 
of the heuristic approach in Section ITl| p^^^ serves as the 



distribution function, and 



(%,v*)V'fe) 



serves as the 



magnitude of tunneling matrix element, and the averag- 
ing of the tunneling matrix element is effectively defined 
such that equation ^ equals equation (10 1. Note that 
the densities of states of both the sample and the tip have 
been reflected in the summation over all states of the tip 
and all states of the sample. 

The above expression for the total tunneling current 
requires diagonalization of the density matrix of the sam- 
ple. To obtain an abstract expression, one needs to 
rewrite it into a form which is a product of a matrix with 
the density matrix of the sample. When only states that 
are electrochemical energy eigen states are used in the 
basis sets, say {|0g)}, with density matrix of the sample 
in the form: 



which is a pure number. As another example, for a Pz 
wave tip state, one has: 



^ip-.Vt) = -—Z—^(^ec,k^^c,p~e.Vt) 



d_ 
dz 



When the pure state Xp is a combination of the s, p, d 
• • • waves, the corresponding operators are additive. 



P0 = 
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this matrix-product form of the expression is: 
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where: 
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Here, the order of stacking (^^;'s for the symbols 0j and ^o 
is the same as that was adopted in writing the density 
matrix p^. Note that the wave functions in the basis 
set {10^:)} are all electrochemical energy eigen functions, 
which makes ^{py^) well defined. 

Although the above results are explicitly written down 



I 

with the basis set {|(/>j:)} which includes only energy eigen 
functions, one does not need this eigen function require- 
ment. The matrices X]p <^o'?^o and J^pfi^c.p — eVt — 

/i)0Q0o have the same transformation rules as other oper- 
ators (and density matrices) in the Hilbert space. Under 
other basis sets, in general one does not have an explicit 
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FIG. 5. Calculation of the tunneling current. The tunneling 
matrix element between one state of the sample and one state 
of the tip is evaluated by a surface integration over the surface 
E between the sample and the tip (see equation (|5|). Refer 
to reference [14] for details. This figure is redrawn in complete 
analogy to Figure 3.1 in this reference. Here the circles rep- 
resent atoms of the sample (left) or the tip (right), and the 
surface E is an imaginary surface in vacuum over which the 
integration in equation ([5| is performed. 



tangled by the (5-symbols, and we were not able to find an 
explicit expression for the measured STP potential fi al- 
though it is implicitly defined by setting equation ( 12 1 to 
be zero. We now introduce some further simplifications 
that make explicit expressions possible. 



D. Explicit expressions for measured potential 

In this subsection we assume that the tip is featureless, 
in the sense that: a) the tip has a uniform density of 
states Nt over the energy range that we are interested in; 
and b) the operator defined in equation (|9| equals: 



^(p.Vt) 



^(^ec,fc'^c,p- 



eVt)J^o 



(14) 



i.e., the p dependence is only in the (5-symbol. 



Under the above simplifications, in equation ( 12 ) the 



sum over p can be explicitly carried out and one obtains: 



Z^ 



Wo = Nt(f)oQ(l)oo 



(15) 



expression that links the matrix elements to the value of 
the sample wave function at a certain point (except for 
the special case discussed shortly in next subsection). 
In equation ( 12 ), the tip and sample properties are en- 



X! fi^c,p - eVt - fi)(fl<Po = Nt$lQ$ooP<p,FD{f^) , (16) 
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where: 
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QfdOj.) = 



l + e'3(^=-A') ' 



(18) 



and Hg is the Hamiltonian of the sample. 

Note that the above equations are basis- 
transformation-independent, hence we do not require 
the basis wave functions to be energy eigen functions for 
the expressions in this subsection. The matrix form of 
the total tunneling current is then: 



/: 
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(19) 



It is noteworthy to compare the above equation to 
equation (Is]). One can immediately see the resemblance 
of the two equations, except that the density of states of 
the sample is implicitly included in the summation of the 
trace operation. One can also see that it is the difference 
of the density matrix that STP measures (with a weight), 
thus in order to express things in terms of effective dis- 
tribution functions, a mean matrix element value should 
be used and that value is expected to change over space 
(as seen in the definition of (/)oo in equation (17l). 

a. Subcase I In this subcase, we consider a situa- 
tion where temperature is high compared to deviations 



from thermal equilibrium due to the applied current. For 
example, in Figure l4] for the heuristic understanding of 
STP, wc need to have fc^T >> eknEo, rather than a 
sharp transition in distribution function implying T « 0. 
We assume the density matrix is very close to the ther- 
mal equilibrium with electrochemical potential fiQ, and 
measured STP potential is fiQ + A/i, then: 



(20) 
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we get: 



2 + e-/3(B-pa) + e/3(^-Mo) ^ ' 
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P tr (t)lo4'OOP'4,.FDitJ'o) 



(21) 



where: 



Q'FD{^J■Q) = 



1 

2 + e-P("--^'o) + e/3(ff=-Mo) 



b. Subcase II In this subcase, we assume that both 
thermal smoothing and nonequihbrium modification of 
the distribution function corresponds to a small recip- 
rocal vector change Afc compared to kp and all other 
characteristic reciprocal vectors. A somewhat cruder ap- 
proximation reads: 



1 



1 



1 + g,3(£;-Mo-AM) 1 ^ pP{E-fio) 




Ar=AK 



Ar=AF, 



hence one obtains: 
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,5(£; - /io)A/i , 



4>m(t>0n {P<t> - P4>,FDilJ'0)) 



where: 



tr </'oo</'ooP"0,-F£)(Mo) 



q" Foim) = S{Hs -Mo) • 



(22) 



The results above constituted a general theory of STP 
that is valid on all length scales and puts coherent and 
incoherent processes on an equal footing through the use 
of the density matrix. The final results from which the 
measured potential is implicitly determined are equations 
^ or (12), which are equivalent. The result is formal, 
however, and not as useful as would be desirable because 
it requires knowledge of the density matrix over the en- 
tire sample - a formidable problem indeed. As we argue 
on physical grounds in the next section, some form of lo- 
cal density matrix would be desirable and seems possible 
to us, although we do not claim to provide a rigorous 
proof. Put more strongly, a formulation in terms of lo- 
cal properties will be necessary in order to make more 
quantitative interpretation of STP result practical. 



IV. LOCAL DENSITY MATRIX 

As argued above, the simple model posed here requires 
including both contacts providing the current and large 
areas of the sample that are far away from the region 
where STP measurement is performed, this is clearly un- 
physical. Moreover, the density matrix of the sample 
is very difficult to calculate numerically due to the large 
number of degrees of freedom. Equally important, ideally 
one would like a theory with which one can uncover rel- 
evant physics of transport without a priori knowledge of 
the density matrix of the sample. A proper theory in this 
case should only involve wave functions in the vicinity of 
the probed area, and involve defects and material char- 
acteristics in roughly the same area. The measurement 



FIG. 6. Schematic of the new basis set that we propose to 
use in order to construct local density matrix of the sample. 
In this figure, only the Af degree of freedom is shown. The 
new basis set is defined with reference to the position of the 
STM tip r, with the states centered around f + Afi. The 
states have finite wave function values at position f when Afi 
is small, and states far away with Afi large have negligible 
wave function values. When the position of the STM tip 
changes, the set {Af} remains the same, hence the central 
locations {f+ Ar^} changes, and the basis set changes. 



result in this case should be determined by these param- 
eters and a proper boundary condition, for example the 
mean current flowing near the region of the sample. We 
argue here that to deal with these issues it is better to 
describe the situation with local density matrices defined 
for each point on the sample, with the entries of the local 
density matrices consisting of only local wave functions. 



A. Outline of local density matrix 

Let us outline here how such a theory might be con- 
structed. To begin, one needs electron wave functions 
that only extend locally. The energy eigen states that we 
used above to calculate the total tunneling current are ex- 
tended throughout the sample, at least in the case where 
there are no defects in the sample. Using this basis set, 
all the states involved need to be counted in expression 
(12), because the wave functions have non-zero magni- 



tude at the center of curvature of the tip tq. In order to 
isolate electrons far away from the measurement area, we 
suggest to go to other basis sets where only local wave 
functions are involved. In principle the density matrix 
of the sample is the same as when we use energy eigen 
states as basis states, however, a new local basis set will 
lead to an approximation that reduces the dimension of 
the density matrix, as we see below. 

We propose the use of a different basis set for each 
point on the sample, see also Figure [6J The basis set 
{l^r+Af fe)} consists of electron wave packets that are 



centered around different positions r + Ari with wave 
vector centered around k. Here, the set {Afi} and {k} 
and the wave function forms are the same for each point, 
but the vector r, which is the position where the STM 
tip probes, changes at different measurement positions. 
The matrices that are transformed from X^p'/'o'^o f^nd 
'^fff{Ec,p — eVt — ^)(t)Q(j)o are thus the same for different 
positions, except for the lateral change in r. While being 
incorrect in a strict quantum mechanical context, with 
this new basis set, one can intuitively talk about electrons 
at certain positions with certain wave numbers. 

Intuitively, the new basis set needs to have Ar,; cover- 
ing the full range of the sample, in order to be a com- 
plete basis set. The matrix form of the density matrix 
of the sample (Cp+AF.,fcl^l?f+Ar,,fe') changes at different 
positions because the basis set changes. 

The simplification from this approach is that density 
matrix elements involving wave packets far away (Af^ 
sufficiently large) should not affect measurement result 
at position r. Note that the density matrix element 
for these states will be finite, but because the region 
of finite amplitude for the wave function is far from 
f, one could expect that both matrices X^n'^o^o s-^id 
Xp/(^c,p ~ eVt — ii)(J)q(J)q in equation (12 1 have zero el- 
ements involving these stateS""^-. With this simplification 
in mind, one might make an approximation and discard 
all states with sufficiently large Ari, and reduce the prob- 
lem to one that only involves electron wave packets near 
to the STM tip position r. This reduced matrix form 
that changes with position is what we define as the local 
density matrix. 

One could thus interpret that what STP probes in the 
relevant non-equilibrium conditions is the local density 
matrix at certain positions on the sample. 

We note that if the total tunneling current was to only 
relate to the local density matrix at one point, the STM 
tip needs to be sharp, with its radius of curvature smaller 
than the lengthscale of variation of the local density ma- 
trix. Only in that case could one choose a surface S in 
equation ([5| that is localized in space. 

The concept of using basis sets with localized wave- 
functions is not new, for example, codes based on linear 
combination of atomic orbital s (LC AO) are common in 
electron transport calculationsSEU, While using these 
basis sets do naturally lead to local density matrices in 
our notion, we would like to point out that these exist- 
ing usages are mostly on small, or even molecular, sam- 
ples; the calculational strategy for a proper calculation 
on large samples, as we will soon discuss, has not been 
worked out. Hence by introducing the concept of local 
density matrix we not only provide a theoretical conve- 
nience, but also pose a new problem in mesoscopic trans- 
port. 



B. Possible further developments 

The local density matrix is defined on a basis set 
{l?f+Af fc)} '^hich consists of only wave packets near 
to the position of STM tip f. The matrices that are 
transformed from Xp'^o'^o and Xp /(-E^c,p— eT4-/i)0j(/io 
in equation ( 12 ) are both determined by transformation 



rules in the Hilbert space, and, once computed, are fixed 
except for the position r for different positions that the 
STM tip probes. In the special case described in subsec- 
tion III D the two matrices above are given explicitly by 



equation (17), no matter what basis set one uses. 



We propose that the approach to calculate the local 
density matrix be different than the path of its defini- 
tion. Since the local density matrix is spatially variant, 
there should be a controlling differential equation that 
determines its variation over space. One can thus pose 
a problem with proper boundary conditions that refiects 
the mean current around the probed area, and need not 
concern the current-providing contacts or the sample area 
far away from the probed area. In other words, although 
the derivation of the controlling equations is based on the 
density matrix of the sample, in the actual calculation 
only the local density matrix is relevant, and one does 
not even need to know the density matrix of the sam- 
ple. The detailed form of the differential equations and 
boundary conditions should depend on what basis set one 
chooses and what further approximations one makes. We 
do not address these computational challenges in this pa- 
per. They require proper theorists. On the other hand, 
below we provide preliminary thoughts of possible further 
simplification procedures. 

The choice of the new basis set is clearly critical in the 
problem. The size of the wave packets in space should 
be small compared to macroscopic lengths in order to 
make the simplification with local density matrix feasi- 
ble; meanwhile this size should be large compared to Xp, 
because of uncertainty principle and the need to keep 
the energy uncertainty small. One can thus imagine two 
limits to simplify the problem: one is to make relatively 
large wave packets, thus including more states in the lo- 
cal density matrix, and making the uncertainty in energy 
small; the other is to make relatively small wave packets, 
which allows to include few states in the local density 
matrix, but at the expense of large energy uncertainty 
for each basis state. 

The matrices that are transformed from Xp'^o'?^o and 

"^ff f {Ec^p — eVt — fi)(j>Q(j)o in equation (12 1 depend on how 
one chooses the new basis set at each point. One can 
calculate these matrices according to the transformation 
rules, or in the case where large wave packets are adopted, 
one may assume that each wave packet is an approximate 
energy eigen state and apply these two matrices in the 
explicit form. 

To simplify the problem without losing important in- 
formation, it is reasonable to assume that all the opera- 
tors J^{f.Vt) are the same for the states in the tip as in the 
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special case of subsection III D viewing the tip as part of 



the instrument and being controlled by the experimen- 
talist. For example, since only nanometer scale signals 
are being probed, one can assume that the tip is an s- 
wave tip with the same constant in the matrix element 
calculation. This approximation was adopted in Chu & 
Sorbello's papep^^. Under this approximation the local 
density matrix approach is also better solidified. 



the sample will become a Fermi-Dirac distribution, i.e., 
suppose the basis set HV'^)} diagonalizes g, one obtains 
from equation (|6|: 



V. LIMITING CASES 

A. Sample in equilibrium 

When the sample does not bear any current, i.e., is in 
equilibrium, after diagonalization the density matrix of 



Ane 



[/(^ec k - M.) - fiEcp - eVt - fit)] 5{E^ 



,fci Ec,p' 



eVt) 



{^ip,Vt)^k) \fo 



(23) 



Note that the Kronecker delta symbol requires the elec- 
trochemical energies for the sample state and the tip 
state to be the same. Under this condition and when 
the temperatures of the tip and the sample are the same, 
if lis > fit, all the terms in the above equation will be 
positive, hence the total tunneling current becomes pos- 
itive; similarly if ^ig < fit, the total tunneling current 
becomes negative. The only possibility for the total tun- 
neling current to be zero is if ^s = Mt, regardless where 
the tip is. This is expected from physical consideration. 
When the temperatures of the tip and the sample are 
different, the above conclusion is no longer true, but this 
is expected as a thermal electric effect. 



B. Homogeneous sample w^ith no defect 

Suppose the sample does not have any defects, namely, 
there is only inelastic scattering and no elastic scattering 
in the sample. Also suppose that the current is uniform 
on the sample. In a jellium model, from translational 
symmetry of the sample (neglecting effects due to the 
macroscopic contacts) , one expects that the local density 
matrix should be the same throughout the probed area, 
except for the effect from a linear electrostatic potential 
distribution reflecting the fact that there is a constant 
current and there exists inelastic scattering. In this case, 
the local density matrix is described by a distribution 
function (no off-diagonal elements), because electrons at 
different eigen states are incoherent to each other. In the 
linear response limit, the distribution function reads: 



N.r ^eE-' 



H' 



3tin.LL/ • I 



(24) 



/(fc,f) ex 



1 



For each direction of momentum there exist an effective 
chemical potential /ig. This distribution function leads to 
the distribution function with electrochemical energy as a 
variant in equation ([2|. In this case, in STP measurement 
the difference of data at different points comes from the 
—E ■ r term. 



C. Landauer resistive dipole 

Suppose there is only one defect on the otherwise ho- 
mogeneous sample, or assume the defects are far away 
from each other compared to Im- Far away from the de- 
fect (s), electrons are not coherent to one another due to 
inelastic scattering events after they are reflected by the 
defect(s). At these positions, the local density matrix 
will again be described by a distribution function. In 
the linear response limit, the distribution function will 
be very similar to equation (24) above, namely: 



f{k,r) ex 



1 



MEs-t^S.p) 



+ 1 



N,r ^ 



-e(/)(r) + ^ + elinVcf) ■ 



(25) 



,/3(Bj-Me,,,) _^ I 



where (pir) is the Landauer resistive dipole potentiaP^l. 
This demonstrates that STP can measure the Landauer 
resistive dipole potential, which is expected from physical 
consideration. Resistive dipo les h ave been observed by 
experiments^. Equations ( 24 ) and ( 25 1 are similar to each 
other, because at the regions where Landauer derived 
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the resistive dipole, the only practical difference is that 
the electrostatic potential is different than the previous 
limiting case due to the existence of the scattering center. 
Here (j){'r) is the dipole field plus the homogeneous electric 
field. 



D. Chu and Sorbello model 

Chu and Sorbello^^ calculated the STP measurement 
result much nearer to a defect than lin. In the paper it 
is assumed that the electrons have a specific distribution 
function in k space before encountering the scattering 
center, and after being deflected by the scatterer, the 
electrons are not deflected again by the same scatterer. 
Interpreting their calculation in our context, they view 
most part of the sample as the current-providing elec- 
trodes that inject the incoming electrons with the spe- 
cific distribution into the "sample" , and the sample size 
is effectively much smaller than /„, centering around the 
defect. Once the electrons enter the sample they remain 
coherent in the presence of elastic scattering. The reser- 
voir case, the background scatterers case and the semi- 
classical barrier case are all describing the distribution of 
the injected electrons. In our context, the density matrix 
of the sample, which is only a very small region near the 
defect, is completely known: 



J2m 



^r 



^1+) 

^k 



(26) 



here 



4->\ 



/ is the purely coherent state resulting from 

one channel of injected electron from the contacts. 

If one considers the case where the part of the sample 
that is not close to the defect is also measured by STP, 
the above density matrix no longer describes the situa- 
tion, that is to say, equation ( 26 ) does not describe all the 
electrons in the system. However, it is conceivable that 
with a local density matrix description, the local density 



matrix very near the defect does look like equation ( 26 1 



Using equation (22) and noticing that f{Ej: — /i^:) 
/(-Bj: - y)~ 5{Ej, - ^J-){^J.^ - m), one obtains: 



SVsTpiro) = - 



Ek\^lt^iroWS{E^, 



m) 



(27) 

which is equation (14) in Chu and Sorbello's paper. 
This illustrates the reason why tunneling matrix elements 
should not be assumed to have the same magnitude in 
writing down the total tunneling current, as the variation 
of this magnitude is the interference term. It can also be 
seen from this expression that the density matrix of the 
sample determines how and where this variation occurs. 



E. Atomic resolution in STP 



As we have emphasized, the tunneling matrix elements 
are dependent on electron wave function values, which 
change over space, leading to the quantum interference 
effects in the end result of STP. In normal STM mode, 
it is also the change of tunneling matrix elements within 
a unit cell that leads to atomic resolution of STM (and 
practically, in order that the atomic corrugations be mea- 
surable, the tip should usually be sharper than an s-wave 
tip, i.e., the operators T acting on the wave functions 
should not be as simple as a pure number, see also dis- 
cussions in Chapter 7 of reference [T4]) . It is then natural 
to ask the following question: does there exist atomic res- 
olution (i.e., corrugations in the end result due to wave 
function value changes within a unit cell) in STP mea- 
surements? 

Part of the answer to the previous question is simple. 
As we discussed in the first limiting case in this section, 
without a current (sample in equilibrium), there shall 
be atomic resolution in STM, but one will see no atomic 
corrugation in STP measurement because it is a constant 
throughout the sample. This demonstrates that if there 
is atomic resolution in STP, it should be different in ori- 
gin than atomic resolution in STM mode. The answer 
to this question when there is a current on the sample 
becomes complicated. To demonstrate the nature of this 
problem, we investigate the STP measurement on a sim- 
ple toy model, namely a one-dimensional tight-binding 
model with only one atomic orbital involved. Readers 
are referred to Chapter 10 of reference [501 for the basic 
setup of tight-binding model. 

We assume that the sample is one dimensional and is 
well described by a tight-binding model with only one 
atomic orbital, the wave function of which is (f>{x). (f>{x) 
is centered around x — 0. The lattice constant of the 
sample is a. In the tight-binding model, the wave func- 
tions of the band is known: 



V'fc(a;)=^e*'=""0(x-na) 



(28) 



where k e (— -, -), which is the first Brillouin zone. For 
our purpose, it suffices to assume that 4'{x) is a s-level 
and to use nearest neighbor approximation, which lead 
to the band structure: 



E{k) = Eq — 7Cos(fca) 



(29) 



where Eq and 7 are two numbers related to the s-level 
energy and AU{x) which is the correction to the atomic 
Hamiltonian. To further simplify things, we assume that 
all operators T are the same for all states (i.e., approxi- 
mation! 14 1 stands, and the expression for total tunneling 



current is equation (19)) for now 
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FIG. 7. Sketch of distribution function with respect to k (up- 
per panel) and E (lower panel), both without current (dashed 
lines) and with a current (solid lines), in the tight-binding 
model without defects under consideration. 



The total tunneling current in this case reads: 
I =^NfY,\F^MxQ)Hl{k)- IFD{k,^i)) , (30) 

k 

where fpoik, fJ-) is the Fermi-Dirac distribution function 
with chemical potential fi (dashed curve in lower panel 
of Figure [7] but has a freedom to move in the horizontal 
direction depending on the value of /x) . The question then 
becomes: does the value of fj. that makes equation (30) 



equal to zero have a spatial dependence (as xq changes 
within the unit cell)? 

a. Xq close to origin When we consider the STP 
measurement results in the vicinity of the atomic orbital, 
only the 7i = component in the definition of ipkix) in 
equation ( 28 ) is relevant, thus the total tunneling current 
becomes: 

/ (X ^ I J-o0(a^o)|' (/(fc) - fFDik, /i)) 



= |J-o(/.(a;o)|'^(/(fc)-/fc(fc,//)) 



(31) 



1. Sample with no defects 

Suppose the sample has no defects, then the den- 
sity matrix of the sample is diagonal under the basis of 
{i'k}- We can use a distribution function /(fc) to de- 
scribe things, and the distribution function with respect 
to energy f{E) can be calculated from f{k). If we ignore 
the thermal broadening, the distribution function with 
respect to both k and energy is sketched in Figure [7] 



Hence there is no spatial dependence in fi. This is in 
contrast to atomicly resolved corrugations in STM mode, 
in which case the change in |Jx)0(2^o)| indicates a spatial 
change in the measured sample height. 

b. Xq dose to a/2 Now consider the case where the 
probed point is not close to the vicinity of the atomic or- 
bital, i.e., at least two (and in nearest neighbor approx- 
imation we consider only two, for example n = and 1) 
components in equation (28) should be considered. The 



total tunneling current then reads: 



/oc 



k 



Ko(; 



{xo) + e''"'Tocl}{xo - a)|' {f{k) - fFD{k,fi)) = {\:Focf,{xo)\^ + \To<l>{xo - a)|') ^ {f{k) - fFD{k,fi)) 



Y, [i^Focj^ixo))* ^Focj^i^o - a)e*'" + :Fo<PM (-^o^l^^o - a))* e-^'^^] {!{k) - fpoiK m)) 



(32) 



r 



It is reasonable to assume that (Jx)0(a^o)) ■^o4'i^o ~ a) 
is real, then the second line of equation (32) equals: 



zero, too. In other words, to this order there is still no 
spatial dependence in the measured STP potential. 



2iTQ(l){xo))* J'Q4){xo-a)Ycos{ka) (/(fc) - fpDik.ij)) ■ 

k 

(33) 
When the current in the sample is small and the system 
is in linear response region, and if /i is the value that 
makes the first line of equation ( 32 ) zero, then the range 
of k for nonzero \f{k) — fpoik, /I)] values consists of two 
very narrow regions centered around —kp and kp, and 
because cos{kpa) = cos(— fci^a), expression (33) equals 



If the current in the sample is large enough, such that 
Afc in Figure [7] is significant, then the above argument 
breaks down and there will be a small corrugation in 
measured STP potential in the atomic resolution. 
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2. Sample with defect(s) 

When the sample does have defect (s), under the basis 
{"i/jfe} the density matrix of the sample has off-diagonal 
elements. As we do not address the calculation of the 
density matrix, we only consider the case in which there 
is one off-diagonal element between states fco and — fco 
(degenerate states and fcp is close to kp). Suppose the 
off-diagonal element is 5 in the (fcoj—fco) entry of the 
density matrix, and 5* in the (— fco,fco) entry. One only 
needs one additional term in equation ( [3l] ) or ( 32 1 in this 
case: 

A/ = To^^koixn) (Toip-kaixo))* 5* + complex conjugate 

(34) 
With the same assumption above for the defect-free case 
that Toipkoixo) (Toip-koixo)) is real, we arrive at: 

A/ =. 2To^Jko (xo) (^oV'-feo (xo))* ReiS) (35) 



a. Xq close to origin In this case, equation (31 1 be- 
comes: 

J2{f{k)~fFD{k,^l))+2Re{5)^0, 

k 

hence there is no spatial dependence in fi. 

b. Xq close to a/2 In this case, the correction to 
equation p2|) includes a term: 



2Re{5) {Fa4>{xt))) To(f>{xQ — a) cos(fcoa) 

X (/(fco)-/FD(fco,/i)) , 

which does lead to a spatial dependence in /i which is 
proportional to S. 

To sum up, when there are static defects in the sam- 
ple, in the transition region between two adjacent atomic 
sites, there will be atomically resolved corrugations due 
to changes in wave function values, but this corrugation 
is a second order effect, as opposed to the first order effect 
in STM mode. If the sample does not have defects, hence 
no interference between states, the atomic corrugation is 
an even higher order effect. 

Though the above results are obtained in a simplified 
toy model, we believe the conclusion bears some univer- 
sality. It also demonstrates the nature of the STP po- 
tential problem, in that to first order the results are of- 
ten trivial, in order to manifest effects from either wave 
function value changes or quantum interference, one of- 
ten needs to go to a higher order. This is also part of the 
reasons that we are not able to provide estimates of the 
size of these effects without calculating the actual density 
matrices. 



VI. SUMMARY AND DISCUSSION 

In this paper, we developed a viewpoint of STP mea- 
surement in the framework of quantum transport. Equa- 
tion ( 12 1 gives an expression of the total tunneling cur- 



the STM tip and the density matrix of the sample. By 
setting this total tunneling current zero, one implicitly 
determines the STP measurement result. The expression 
is in a matrix product form such that a basis-set free def- 
inition can be made. With some approximations, we also 
obtained explicit expressions for measured STP potential 
in the featureless tip case. 

In order to get rid of the unphysical requirement of 
material characteristics far away from the STP probing 
area in writing down the density matrix of the sample, we 
proposed to use local density matrices, defined in Section 
IV as a varying parameter over space, which is controlled 
by a certain differential equation and a set of boundary 
conditions. This allows relating STP to sample proper- 
ties near the region where STP is performed. 

In Section |V] we provided limiting case calculations 
that demonstrate the usage of the theoretical descrip- 
tion, mostly with given or assumed density matrices of 
the sample. In particular, we found that with the toy 
model of one dimensional tight-binding metal, atomic 
resolution in STP measurement is a high order effect. 
Following the derivation of this particular example, one 
can find that the null results are partly due to assumed 
symmetry in the Fermi surface (e.g.. Equation (33) being 
zero under the adopted assumptions) or in the operators 
T acting on the sample wave functions. We speculate 
that asymmetry in the Fermi surface with respect to fc, 
or in the tip wave function (hence operator J^) with re- 
spect to energy, could lead to lower order effects in STP 
measurement. 



Comparing equation ( 12 ) with, for example, equation 



rent which is dependent on the distribution function of 



6.6) in reference fTSl one can immediately see the anal- 
ogy. The matrices Y.p$o$o and Y.pf{Ec,p-eVt-n)$l$o 
are in the same place where the scattering matrices are. 
Although one needs to keep the subtle difference between 
correlation function and density matrix in mind, these 
two matrices can be viewed as a special kind of scatter- 
ing matrix. 

The STM tip as a special contact provides scattering 
which is described by a scattering matrix related to the 
matrices mentioned above. This implies that it will also 
act as an electron source, the strength of which is depen- 
dent on the strength of scattering. The electron source 
effect will change the density matrix of the sample, or 
the local density matrix, as the tip moves. However, in 
the limiting case that the strength of scattering is small, 
i.e., the probe is weakly interacting with the sample, STP 
measurement result is not expected to change much. STP 
does have this advantage of perturbing the sample mini- 
mally, compared to other contact-based probes, for exam- 
ple, using conducting cantilevers as the fourth electrode 
and using an atomic force microscope to scan on the sam- 
ple. Conceptually, a good criteria for minimal source ef- 
fect would be that the tip should be far away from the 
sample that it is in the far tail of the exponential decay of 
electron wave function for sample states, though this will 
correspond to a large tunneling resistance and the reso- 
lution of the STP measurement in this case is expected 
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to be low due to large Johnson noise. Nevertheless, it 
is a valid starting point to work under this assumption 
to simplify the problem. When the sample is small com- 
pared to its characteristic transport lengths, the effect 
of the tip is generally large and these effects are alrea- 
day explicitly included in the formalism of works on this 
cas^ 
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